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We consider the superconducting proximity effect in a helical Luttinger liquid at the edge of a 2D 
topological insulator, and derive the low-energy Hamiltonian for an edge state tunnel-coupled to a 
s-wave superconductor. In addition to correlations between the left and right moving modes, the 
coupling can induce them inside a single mode, as the spin axis of the edge modes is not necessarily 
constant. This can be induced controllably in HgTe/CdTe quantum wells via the Rashba spin-orbit 
coupling, and is a consequence of the 2D nature of the edge state wave function. The distinction of 
these two features in the proximity effect is also vital for the use of such helical modes in order to split 
Cooper-pairs. We discuss the consequent transport signatures, and point out a long-ranged feature 
in a dc conductance measurement that can be used to distinguish the two types of correlations 
present and to determine the magnitude of the Rashba interaction. 

PACS numbers: 74.45. +c, 71.10.Pm, 73.23.-b 



I. INTRODUCTION 

The helical edge states of a 2D topological insulator 
(TI) consist of a Kramers pair of right- and left-moving 
electron modes of opposite spin situated inside the bulk 
gapPS and they have so fa r bee n observed in HgTe/CdTe 
quantum wells (HgTe-QW) 2 * 5 * 6 l In 3D topological insula- 
tors, the edge states cover the surface of the material and 
consist of a single- valley Dirac cone with spin-momentum 
locking, which leads to unique electromagnetic proper- 
ties and quantum interference effects 7 . In both 2D-TI 
and 3D-TI the coupling of spin and orbital motion can 
lead to interesting effects when combined with supercon- 
ductivity. Superconducting correlations induced by the 
proximity of a singlet s-wave superconductor can inside 
the TI obtain a p-wave character, which can be used to 
engineer Majorana bound statesP*^ A somewhat sim- 
ilar induction of non-conventional correlations has also 
been proposed to occur in other semiconductor systems 
in the combined presence of the spin-orbit interaction and 
superconductivity.^ 

When the edge state of a 2D-TI is coupled to a sin- 
glet superconductor, the transfer of electrons between 
the systems can, first of all, induce singlet-type prox- 
imity correlations between electrons in the right and left 
moving modes (the H — channel) W This already leads to 
several effects of interest. For instance, the helicity of the 
electron liquid lifts the spin degeneracy and enables Ma- 
jorana states,^ causes Cooper pairs to splitf^and affects 
transport properties.^ Tight-binding calculations study- 
ing the pair amplitude have also been made^. There is, 
however, also a possibility of inducing correlations only 
within the right-moving (or the left-moving) channel at 



a nonzero total momentum (the ++ and channels). 

Such a channel is not forbidden by symmetries in the 
problem: due to the spin-orbit coupling, the spin axis 
of the edge state is not necessarily constant, so that the 
electrons forming a Cooper pair singlet can both enter 
the same mode on the TI edge, even when spin is con- 
served in the tunneling process and time-reversal symme- 
try is present. In HgTe-QW, a non-constant spin axis can 
be induced externally by the Rashba spin-orbit coupling 
that breaks inversion symmetryP Momentum conserva- 
tion is required to be broken, but this can occur e.g. due 
to inhomogeneity or a finite size of a tunneling contact. 
Moreover, unlike in metals, in 2D-TI the momentum non- 
conservation can in principle be made arbitrarily small by 
tuning the Fermi level near the Dirac point (k = 0). 

A straightforward way to probe the existence of su- 
perconducting correlations is to observe the Josephson 
effect or other interference effects that can be mod- 
ulated with superconducting phase differences. The 
Josephson effect has been studied previousl y i n vari- 
ous one-dimensional Luttinger liquid systemsP^HII] The 
finite-momentum channel has, however, received limited 
attentionpSO and is usually negligible. As shown below, 
certain experiments with superconducting contacts at- 
tached to the helical edge states can nevertheless probe 
such microscopic aspects of the tunneling, including the 
role of the Rashba interaction. 

Here, we first derive a low-energy Hamiltonian de- 
scribing the superconducting proximity effect in the edge 
states of a 2D TI coupled to a conventional superconduc- 
tor by tunnel contacts. We use it to find the signatures 
of both types of tunneling events in a transport exper- 
iment. Because of the reduced number of propagating 
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modes in the helical liquid, correlations within the same 
channel occur at a finite momentum and, as in chiral 
liquids,^ are affected by the exclusion principle. It turns 
out that although this component of the proximity effect 
gives a negligible correction to the dc Josephson effect, 
in the NS tunneling conductance [see Fig. [ljc)] it man- 
ifests as a long-ranged interference effect, oscillating as 
a function of the superconducting phase difference, and 
unlike the H — part, is not exponentially suppressed at 
length scales longer than the thermal wavelength. The 
ratio of the contributions of the two possible channels 
scales as SG++/SG+- oc (z /hv F ) 2 (k B T/M) 2 e 2 ^ hVF 
(in the noninteracting case), where zq characterizes the 
strength of Rashba interaction, vp is the Fermi velocity 
of the edge channels, M is the energy gap of the TI, and d 
the distance between two superconducting contacts form- 
ing the interferometry setup. The amplitude of the effect 
is proportional to the amount of spin rotation achieved 
by Rashba interaction, and the quadratic temperature 
dependence is due to the exclusion principle. We also 
discuss how e-e interactions modify this result. 

This paper is organized as follows. In Section[llJ we in- 
troduce the model for the helical Luttinger liquid (HLL) , 
the coupling to the superconductors, and the electronic 
structure of HgTe-QWs. Section [III] discusses the effec- 
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tive low-energy Hamiltonian, and Section IV transport 
signatures in the dc and ac Josephson effects and the NS 
conductance. Section [V] concludes the manuscript with 
a discussion on the results and remarks on experimental 
rcalizability. 
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FIG. 1. (a) The setup considered: a 2D topological insulator, 
whose edge state is coupled to two superconductors via tun- 
nel contacts. In response to phase ip or voltage V differences 
between the superconductors, Josephson currents can flow via 
the edge channel, or one can probe the NS transmission by in- 
jecting current from the superconductors to the edge channel. 
If the spin axis of the edge state is not constant spatially and 
as a function of energy, a Cooper pair singlet can enter the 
edge state in two possible ways: either electrons enter modes 
propagating to opposite directions (left) or the same direc- 
tion (right). In the latter case, exclusion principle requires 
temporal (or spatial) separation of the two, which can be of 
the order of the superconducting coherence length £ = Kvf / A 
still preserving the correlation, (b) Configuration for the mea- 
surement of the Josephson current, (c) Configuration for the 
measurement of interference in the NS conductance. 



II. MODEL 

We consider the setup depicted in Fig. [I] The edge 
states of a 2D-TI are coupled to two superconducting 
terminals via two tunnel junctions. Below, we in general 
assume that the distance d between the contacts is longer 
than the superconducting coherence length £. 

The left- and right moving edge states |+,x) and 
|— ,x) have a linear dispersion, and are described by the 
bosonized HamiltoniarP 



Ho 



&xu\ 9 - 1 {d x $y +g{d x cf)f 



(1) 



where the Fermi field operator 
ip a (x) (2Tra )- 1 / 2 U a e lak " x e t ' l 'o^ 
(2 7 Ta )- 1/2 U a e la ^ FX+ ^ i> ^ +l ^ T ^ x \ the stan- 
dard boson fields $(x), cp(x) satisfy [0(ac), ^{x')] = 
(i/2)sgn(ir. — x'), and U± are the Klein factors. 
u = vp/g is the renormalized Fermi velocity. Here and 
below, we let h = ks = e = 1, unless otherwise men- 
tioned. The parameter eto is the short-distance cutoff. 
In the noninteracting case, the Luttinger interaction 
parameter g = 1, and with repulsive electron-electron 
interactions one has g < 1. 

The coupling to the superconductors is modeled with 



a tunneling Hamiltonian 



a=±,c'=t4 



dx dV t au > (x, (x)ip s *> (7) + h.c. . 

(2) 



where the tunneling amplitude t aa i(x, f 1 ) describes the 
tunneling from the state \o' , r") in the superconductor to 
state | a, x) in the edge mode. For what follows, it is use- 
ful to introduce also the corresponding one-particle oper- 
ator Iit, in terms of which, t arT >(x, r 1 ) = (a, x\hr\cr' , f) . 
The momentum k along the edge is a good quantum num- 
ber for straight TI edges, and we define the state la, x) in 
is the momentum representation: \a,x) = e~ lkx \a, k), 



where \a, k) is the edge eigenstate with momentum k and 
propagation direction a = ±. 

We assume that the Hamiltonian is time-reversal sym- 
metric, which implies that the tunneling operator in gen- 
eral satisfies T/i^T -1 = hx- Here, we choose the phases 
of the wave functions so that the time reversal operations 
read %\a' ,?) = a'\-a' \r>) and %\a,k) = a\-a,-k). We 
also assume that the tunneling is spin-conserving, that 
is, written in terms of real electron spin states in the TI 
and the superconductor, we have (a, r\hr\~ &i O = 0. 

To describe tunneling to HgTe-QWs, we need some 
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knowledge of the structure of the edge states. This can 
be obtained from the four-band model used in Ref. [2] 
In this approach, the low-energy properties of the TI 
are described using a 2D envelope function in the basis 
of four states {\E1+), \H \E1-), \H1-)} localized in 
the quantum we ll. 121221 The edge states at the boundaries 
of the TI can be solved within this four-band model J^S for 
which we give a full analytical solution in Appendix [A"| 

We assume the terminals are conventional spin-singlet 
superconductors. As usualpl they are characterized 
by the correlation function F(fi, a±, n; f 2l o 2 , t 2 ) = 
(T[ip (Tl (f*i, Ti)ip a2 (f 2 , t 2 )]}q that has a singlet symmetry 
F(fi,a ll Ti;f 2 ,a 2l T 2 ) = aiS ai -^i^fi, n; f 2 , r 2 ). In the 
bulk, the correlation function obtains its equilibrium 
BCS form, which in imaginary time can be written as 



F(ri,r 2 \w) 



d 3 k 
(2tt) 3 



a — i(f-L—T t 2)-k_ 



A 



^ + & + \M 2 



(3) 



with £k = k 2 /(2m) — [i the dispersion relation and A the 
gap of the superconductor. 



III. EFFECTIVE HAMILTONIAN 

Integrating out the superconductors using perturba- 
tive renormalization group theory (RG) and consider- 
ing only energies \E\ -C |A| reduces the Hamiltonian 
Hq + Ht of the total system to one concerning only the 
one-dimensional edge states: 



H = H a + I dx T+_ (x)fl>+ (sty- (x) (4) 

+ ^++( x )^+( x )^+( x + a) 

+ T (x)ip-(x)i/)-(x + a) + h.c] . 

Here, T a p describe the coupling to the superconductor, 
and a = hvp/A is the new short-distance cutoff in the 
theory. Details of the derivation are discussed in Ap- 
pendix [B] 

The coupling factors in the noninteracting case (g = 1) 
are given by the expressions (see Appendix [B] for general 
discussion) : 

T ++ (x) = U dMdM 5>^F(r'i,^;0) (5) 



xAvp 1 d k P ++ (j + k,r / 1 ; J - -/.-.,■' ' 



K 

K 
~2 



2 ) \k=0 : 



where F is given in Eq. ([3|, and 

T + _(x) = n f d 3 r[d 3 r> 2 ^ e iK 'F(^, r' 2 ; 0) (6) 

J K 

K K 

x^+-(y ~ k Fl r[; — + k F ,r' 2 )* , 

with r = in the presence of time reversal sym- 
metry. The main contributions should arise around 




FIG. 2. Schematic depiction of the edge mode spectrum 
and the states participating in 2-particle tunneling from the 
Fermi level in the superconductor to states in the TI. (a) In 
the -I — channel, the momenta corresponding to energy e are 
fci = kp — e/hvF, k2 — — A;_f — e/hvF- (b) In the ++ channel, 
one has k\ — + e/hv F , k% = —k F — e/hv F . Singlet pair 
tunneling into the ++ channel can occur if the electron spin 
axis is different for the states at ki and ki. 



K = —2k F for r ++ (due to the 2kp oscillations in the 

Fermi operators), and around K = for 

The coupling is proportional to the factor 



P ai a 2 (h,r[;k2,r 2 ) = [t aii (k 1 ,r' 1 )t a ^(k 2 ,f' 2 ) 
- t ai ^(k 1 ,r' 1 )t a2i (k 2 ,r' 2 )\ ■ 



(7) 



which describes two-particle tunneling of a singlet, from 
two points and rJ> in the superconductor, to momen- 
tum states |ai,— fci), \a 2 ,—k 2 ) in the TI edge modes 
(cf. Fig. g. Here, W^f) = / dx e lkx t aa , (x, f) = 
(a, — k\fvr\a' , f) is a Fourier transform of the tunneling 
matrix element. 

One can also verify that in the absence of interactions, 

the expression for the amplitude coincides with the 

leading term in the zero-bias conductance in the normal 
state, up to a replacement F m- 2(nvp)~ 1 lmG R . Within 
a quasiclassical approximation in the superconductor,^ 
one then hnds a relation to the normal-state conductance 
per unit length, g(x), of the tunnel interface: 



"+-{x) ~ ^hv F R K g(x) 



hv F Rk 

~h AR 



N 



(8) 



Such a relation is typical for NS systems, and connects 

the amplitude r_| to observable quantities. The latter 

expression assumes the total resistance Rn is uniformly 
distributed in a junction of length l F . When the inter- 
face resistance decreases, the effective pairing amplitude 
r grows — and although not included in our perturba- 
tive calculation, one expects that this increase is cut off 
when the effective gap reaches the bulk gap of the super- 
conductor, = A. 

Unlike , the amplitudes do not have a 

direct relation to the normal-state conductance, and they 

depend on the factors P++/ , which are proportional to 

the spin rotation between the states involved in the pair 
tunneling (see Fig. [2]). Estimating this factor is necessary 
for determining how large the same-mode tunneling is in 
a given system. 
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A. Two-particle tunneling 

Making use of the time-reversal symmetry, it is possible 
to rewrite the P factors in a more transparent form: 



P aia2 (ki,r[;k2,r' 2 ) = (a 1 ,-ki\Z(r\,r 2 )%\a 2 ,-k 2 ) , 
Z{^ 2 )^h T [l^{\^){^ 



"' IA|<!2. -'.'2/ 

(9) 

f 2 )(^\)}h T , (10) 



where \ a is the identity matrix in the spin space of the 
superconductor. Unlike the starting point, this expres- 
sion is explicitly independent of the choice of the spin 
quantization axis. We also note the symmetry: 

P axa2 (ki,r[;k 2) r 2 ) = -P a2ai (k 2 ,r' 1 ;k 1 ,r 2 ) , (11) 

following from the definition Eq. 

We can now make some remarks on the possibility 
of ++ tunneling. First, suppose that the state \a, k) 
describes an electron wave function with a fixed fc- 
independent and spatially constant spin part, and that 
the tunneling is spin conserving. In this case it is easy 
to see that P ++ = 0, as the inner product of a spinor 
and its time reversed counterpart vanishes. Such a situa- 
tion is realized, for instance, within the plain Kane-Mele 
model: 1 Breaking such conditions can, however, lead to 
P ++ 7^ 0. We demonstrate in the next section that this 
can occur in HgTe-QW. 



B. Effect of Rashba interaction in HgTe/CdTe 
quantum wells 

We now discuss a simple model for tunneling into the 
helical edge states of a HgTe-QW, taking spin axis ro- 
tation from the Rashba interaction into account. We 
make the following assumptions: the tunneling is spin- 
conserving and local {(r\hT\r') oc S(r — f)] on the length 
scales of the four-band model. This results to all contri- 
butions to P ++ coming solely from the Rashba mixing. 

While we cannot estimate the actual values of P_| or 

P++ within this simplified a model, we can study their 
relative magnitudes, which is now determined by the low- 
energy four-band physics only. 

Under the locality and spin-conservation assump- 
tions, the tunnel matrix element introduced above ob- 
tains the following form in terms of envelope spinor 
wave functions W in the four-band basis {\j}} = 

{\Ei+),\Hi+),\Ei-),\m-)y. 

P aia2 (k 1 y i ;k 2 ,f' 2 ) = * 01 - kl (x'ny'xy Z^iA) 

x%i> a2 ,-k 2 {x' 2 ,y' 2 ) + [r< x ^r< 2 ], 

(12) 

Z&rthj' = 01Ml«r ® K)(r> 2 \]h T \f) . (13) 

Time reversal for the four-band spinor reads T = ~iT y K 
with K the complex conjugation, and the r matrix acts 
on the Kramers blocks (+, — ). For simplicity, we use now 
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FIG. 3. Projections z(k) of the additional spin-orbit cou- 
plings Rq and To on the edge state basis, for the parameters 
of Ref. El with M = -lOmeV. The dashed lines indicate lin- 
ear approximations Z\ ~ Zo,ife and zz ~ zo,3& with 2o,l = 0.03 
and zo,3 = 0.03nm -2 . 



a length-scale separation between the scales appearing in 
the four-band model and the atomic ones (tunneling 
Z, kp^s m the superconductor, unit cell). We consider 
only the long-wavelength part of P, and replace Z with a 
constant describing the tunnel coupling to the quantum 
well basis states, obtained by averaging it together with 



F [cf. Eqs. (J5j>, ^} over r[ and 



[Z{^ 2 ) + Z{r> 2 , ^)}F{^ 2 . 




o 



23ft) 
-Vft) 



(14) 



with A and B real-valued. This form follows from the 
time reversal symmetry and hermiticity of the matrix el- 
ements of the operator in Eq. ( 10 ). We have also assumed 



here that the decay length for the F function (~ kp S ) 
is short on the scales of the 4-band model. Finally, we 
for simplicity neglect the coupling to the HI band, and 
set B = C = 0. For a lateral contact (SC on top of 
HgTe-QW), the main tunnel coupling is expected to in- 
volve the El band, which extends deeper^ into the CdTe 
barrier than HI. Including additional couplings would 
however cause no essential qualitative differences in the 
estimated ratio between the ++ and +— terms. 

Without additional spin axis rotation from the Rashba 
interaction, the edge states are in separate Kramers 



blocks (see Appendix [AJ), * + oc ($+,0) and ^_ oc 

(0, $-), and we can see that P ++ = whereas P^ oc A. 

Note that a contribution proportional T> does not arise: 
the unperturbed edge state wave functions are both pro- 
portional to the same constant real-valued spinor, &± oc 
X, so that the Independent contribution would be pro- 
portional to X i a yX — 0- This structure also implies that 
contributions proportional to T> do not arise in the lead- 
ing order of the Rashba coupling. 

Rashba and other related spin-orbit interactions in the 
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four-band model can be represented as 



H R = 











h R = 



iSnk 2 _ 



—Rf)k_ <5 "1" 

-S - iS Q k 2 _ Tokl 



(15) 



where k± = k ± ik y . For the QW parameters used in 
Ref.EU R « -15.6nm 2 x eE z , iS w -2.10nm 3 x eE z , 
and T w — 8.91 nm 4 x eE z , where E z is the electric field 
perpendicular to the QW plane. The model could also 
include the bulk inversion asymmetry terms <5P^ 

To obtain the effect of the Rashba interaction on the 
wave functions, we find the low-energy eigenstates of H = 
Ho+Hr numerically. For given k, this is a 1-D eigenvalue 
problem in the y-direction, which can be discretized and 
solved by standard approaches. Analytical results can be 
obtained by perturbation theory in Hr restricted to the 
low-energy subspace spanned by the unperturbed edge 
states. For typical experimental parameters, the whole 
wave functions ^ however turn out to have a significant 
component also in the continuum of bulk modes above 
the gap, which is not adequately captured by such an 
approach. Our estimates for the matrix elements P a p 
below are therefore based on the numerical solutions for 
the eigenstates. 

However, qualitative understanding can be obtained on 
the basis of the model restricted to the low-energy sub- 
space. Projecting Hr to this basis (see Appendix |A|, we 
find the effective low-energy Hamiltonian of the system^ 



i , -i{R Wi{k) + T w 3 (k)' 
Hr - 1 c.c. 

w 1 (k)=xf I dyf +tk (y)[k + d y ]f- >k (y) 

-OO 



(16) 
(17) 
(18) 



where $±.fc(y) = f±,k(y)x- This result is valid to the 
leading order in tiR. The constant and quadratic in k 
terms (proportional to 6 and So) give no contribution, as 



X^yX 



0. Using typical HgTe-QW parameters,^ the 
integrals evaluate to u>i(fc) ~ Wo,ik and w^k) WQ^k 
near the Dirac point, as illustrated in Fig. [3} The prefac- 
tor u>o,i ~ 0.03 is essentially independent of the mass pa- 
rameter M, and w ,3 ~ 0.03nm~ 2 x (|M|/10meV). Note 
here that the matrix element 0.03i?o& of the Rashba in- 
teraction with the edge states is significantly smaller than 
the Rok± appearing in the bulk Hamiltonian. The 2x2 
effective Hamiltonian yields the wave functions: 



+,k 



2v F ^~ 




(19) 



where wq — RqWq,i + ?oWo.3- The Rashba interac- 
tion mixes the two Kramers blocks, but in the lead- 
ing order does not modify the energy dispersion. Al- 
though the mixing angle of the $±,fc spinors is indepen- 
dent of k, the total four-band spinor is not: the decay 
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FIG. 4. The El- component of |*+,-k| at k = 

\M\/{61wf). It is linear in the Rashba parameter wo/vf, 
provided \wo/vf\ < 1- Shown are numerical results (solid 
line), the component in the unperturbed edge state subspace 
(dashed line), and the result of Eq. |T9| l (dotted). The be- 
havior for y > 20 nm depends mainly on the linear in k 
Rashba term Ro- Here, M = — lOmeV. Inset: results for 
M = -50meV. 



lengths l/Ai/ 2 (fc, a) of &±,t in the y-direction depend on 
k and are different for the a = + and a = — states: 
time-reversal symmetry only guarantees Xi/2(k,+) = 
Ai/2(— k,— ). This makes the electron spin axis to ro- 
tate both spatially and with energy e, which ultimately 
is required for a finite P++. 

For comparison, we show in Fig. [4] the El— compo- 
nent of the numerically computed total edge state wave 
function ^+.k{y), and its projection to the low-energy 
subspace, which can be seen to match Eqs. ( 19 1 to a 
very good accuracy. The El— component is propor- 



tional to the Rashba coupling and contributes to P++. 
As is clearly visible in the figure, neglecting the bulk 
states underestimates the total amount of spin rotation, 
for experimentally relevant parameters. For a larger (but 
unphysical) value for the gap \M\, the low-energy theory 
works slightly better, as visible in the inset of Fig. [4] 

We can now estimate the relative order of magnitude 
between P++ and P^ within this model. From the re- 
sults above, one can see that the representative quantities 
to be compared are 



and 



ihvF k ++ ^ 2 



? 2 ) = P + -{*-k F , 



K -K 

2 



(20) 



-f K 7 



(21) 



In Fig. [5] we show the ratio of these amplitudes for f\ — 
7% = (0, y) (i.e., the value at a distance y from the edge). 
The C++ amplitude increases when the energies of the 
edge states involved approach the TI energy gap edge. 

The general order of magnitude of the c ++ u factors 

can be estimated to be of the order 



C-H 



'-■() 



\M\ hv F 



(22) 
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FIG. 5. Relative magnitude of the two types of tunneling as 
a function of location and parameters, at K = —2k F for P++ 

and K — for , as obtained from numerically computed 

ty±,k- The scaling with the electric field is given by zq = 
leVnm x _B z /(500 mV / nm ). Solid and dashed lines indicate 
hv F k F = ±\M\/(6hv F ), i.e., energies E = E oi[!LC q= \M\/6 ^ 
(0.75=F0.15)|M| for the states involved in the c+ +/ /+_ factors. 
For Zo/vf <■ 1, C++ is linear in zq. 

A similar relation is then expected also between the T++ 

and r_j factors for surface contacts to area near y — 

0. Here and below, we characterize the strength of the 
Rasbha interaction with the quantity zq = 1 eV nm x 

nm I . 

Using the above results for the order of magnitude of 

C++ and C-| we find (see Appendix B) the estimates for 

the general case with e-e interactions: 

r + _ ~ (aoAJ^-T'l 1 (23a) 
r ++ ,i (O0 A>^--^n-, (23b) 

corresponding to cutoff a = Hvf/A. The relation Q 
for the noninteracting value r^" 1 ~ (Hvf/It)Rk /(4-Rjv) 
fixes the magnitudes relative to experimental parameters. 

With finite electron-electron interactions (g ^ 1) in the 
helical liquid, all the effective tunnel rates obtain identi- 
cal scaling in the original short-distance cutoff ag. This 
reflects the renormalization of the single-particle tunnel- 
ing elements t aa by the electron-electron interactions. 

We can also estimate the Rashba coupling factor ap- 
pearing in T++. With a typical TI gap \M\ ~ lOmeV ~ 
100 K, we see that the factor of A/|M| can be made of 
the order of 0.1... 0.2 with conventional superconduc- 
tors, and can be even larger for smaller TI gaps. The 
second factor is zq/Kvf ~ i? 2 /(100 mV /nm), and as vis- 
ible in Eq. ( [19] ), measures the rotation of the spin axis 
caused by the Rashba mixing. An upper limit for the 
field that can be applied in practice is likely of the or- 
der E z ~ 100 mV /nm, as for fields larger than that, the 
potential difference across the QW becomes compara- 
ble to the energy gap of the barrier material (CdTe). 
Based on this we find an estimate for the achievable ra- 
tio, r++ ~o.i...ir+_. 



Finally, let us remark that tunneling that is local 
in real space, (flftrK) ^ 8{f — f), does not lead to 
tunneling that is local in the edge state Hamiltonian, 
t aa (x,r l ) cx S(x ~ x'). This follows in a straightforward 
way from the extended 2-D nature of the edge states 
and the k x , k y mixing due to the spin-orbit interactions: 
(a,x\h T \a,r) fx £ fe e^" 35 '^^^, y', z')*. If the spa- 
tial profile y of the wave function has /c-dependence on 
the scale k , the sum resembles a rounded 5 function of 
width fc^ 1 . For HgTe QW edge states, k^ 1 - hv F /\M\ 
is a low-energy length scale. Because of this, a pointlike 
contact to a superconductor can produce a finite P++, 
even though assuming t(x, f") oc 5(x — x') in Eq. |7|) leads 
to the opposite conclusion. 



IV. TRANSPORT SIGNATURES 

To study the experimental signatures implied by the 
above model, we consider the transport problem in the 
setups depicted in Fig. [T] There, two superconducting 
contacts are coupled to a helical liquid, whose potential 
is tuned by additional terminals at the ends. There are 
three related transport effects one can study here: the 
equilibrium dc Josephson effect, the ac Josephson effect, 
and the NS conductance. 

We consider a general nonequilibrium case of a time- 
dependent pair potential A(i) = |A|e l<Pl ^^ in the left 
contact and A(t) = lAle 1 ^ 2 ^ in the right one, with 
ip^t) = sp /2 + 2V 1 t and <p 2 (t) = -<p /2 + 2V 2 t. In 
Eq. Q, the factors T inherit this time dependence. We 
also assume that only sub-gap energies are involved in 
the transport, so that the quasiparticle current to the 
superconductors remains exponentially suppressed by the 
superconducting gap. 

The current is obtained as an expectation value of a 
current operator / = i[H,N] where N is the particle 
number in the HLL. From the effective Hamiltonian, we 
identify 




where Isi and 1 32 must be interpreted as the parts cor- 
responding to currents injected through the interfaces at 
SI and S^. The sums over af3 run over H — h, H — , and 

Considering only the Cooperon terms [cf. Fig. [7Fa)], 
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using perturbation theory up to second order in T we find 
IjM*) = -SlmY" / da* / daaX^Cn.xa) (27) 

a p J SI JS2 
poo 

x e^'W / di' e -^(*-*') Imlxa^X! - x 2l *')] , 
Jo 

Ij,S2(t) = Ij,si(t)\ipi^ipi , (28) 



where 



x a0 ( Xl) x 2 ) = e 2lk ^ a+ ^- x ^r aP ( Xl )r a0 (x 2 y (29) 

/ e i0«(a;,t) e i0u(x,t) e -i^ Q ,(O,O) e -i^(O,O)\ 



Xa^(a;;*) 



(27ra) 2 



(30) 



The a/3 = H — component of the current coincides with 
the result obtained in Ref. [T5] Note that the terms in- 
cluded here contain the leading order of the dependence 
in the phase difference ipi — ip 2 . 

The above correlation functions can be evaluated via 
standard bosonization techniques!^ 

X * a (x,t) = (27ra)- 2 B a (x 7 ty+^ 1+2 B_ a (x,ty + 3~ 1 - 2 , 

(31a) 

X+-(x,t) = (2na)- 2 B + (x,t) 1 ^B_(x,t) 1 ^ , (31b) 

B± (x, t) = . ~ mZ . , (31c) 

sinn[z(w — 1a =F x)\ 

where z = irT/u. 

In the noninteracting case (g — 1), we can evaluate the 
time integrals analytically, to order 0(a 3 ): 

Ij,si(t)= [ tei I dx 2 [j+++^ 1 +j+5 1 ] (32) 

J SI JS2 

3j,Si = l ~^V 2 [(V 2 /A) 2 + 4n 2 (T/A) 2 } (33) 

x cosT^o + "2Vit - 2V 2 (t - \xi - x 2 \/v F ) + 4>oj , 
jj,si = 0, (34) 

- l ^." z n (35) 

x sin^o + ZVyt - 2V 2 (t - \xi - x 2 \/v F ) \ , 

where (f>o(xi,x 2 ) = 4fc_ F (xi-a;2)+arg[r ++ (a;i)r ++ (.'E2)*] 
is a dynamical phase shift. 

Below, we discuss the implications of these results first 
at equilibrium and then at finite biases. 



3j,si 



A. Equilibrium 

At equilibrium, the leading contribution to the super- 
current comes from the H — channel. As shown in Fig. |6j 
the supercurrent is finite at zero temperature, and de- 
cays exponentially as the temperature is increased above 




1 2 3 
nTdlhvf 



FIG. 6. Equilibrium Josephson current between the su- 
perconducting contacts, relative to its noninteracting zero- 
temperature value 7o. (a) Dependence oc A 2-2 ^ 9 (T/A — 0.1) 
on the interaction parameter at T -C fof / d. (b) Temperature 
dependence for different values of g. 



hvF/d, in a way that depends on the strength of electron- 
electron interactions. The qualitative features are the 
same as those found in Ref. IT51 

The contribution from the ++ and channels to the 

equilibrium current is not more relevant than H — even in 
the interacting case, unlike in Ref. Q2B Based on scaling 
dimensions in the effective Hamiltonian (dim?/>+'0- = 
g , dimi/j + t/j + = g + g^ 1 ), one finds the scaling 1^ oc 
(£/A) 2 /9- 2 and J++/~ oc {E / Af^+a' 1 )-^ for the low- 
energy scale E = max(T, Vp/d), which implies that 1^ 
will be more relevant than whatever the inter- 
action parameter. This difference arises from the ex- 
clusion principle, which makes the + + / channel 

less favorable for the supercurrent, although note that 
with decreasing g (larger repulsive e-e interaction), the 

(+ + / ) contribution grows relative to the H — one. 

However, as noted in Section|Bj the scaling with the bare 
short-distance cutoff a as opposed to A -1 is identical for 
r _|__|_/ and 



B. Nonequilibrium 



When the superconductors are biased with a finite volt- 
age, currents generically start to flow between all the ter- 
minals, and they may also be time dependent due to the 
ac Josephson effect. To fully understand these effects, it 
is illuminating to compute the spatial distribution of the 
currents in the system. 

The spatial dependence of the currents in the helical 
liquid can be obtained by making use of the following 
expression for the current operator / = ^r=d x (/)(x,t) in 



the Heisenberg picture (see App. [c])!^ 



I(x, t) = Iq(x, t) + v F 
1+5 



Ax 1 dfr 



o>=± 



(36) 



D (^x ^ t ^ x j ^ ) — q/ ^«z> 3> j i ) 



j a {x',t') = 



25 
5(j> a {x') ' 



2.9 



(37) 



This applies to any Hamiltonian of the form if = ifo + 
^(t), where -f/o is the bosonized Hamiltonian in Eq. ([I]); 
Jo is the current operator evolving in time with the unper- 
turbed Hamiltonian H . The operator j a (x' ,t') can be 
interpreted as the current density injected to the mode 
a = ± at position x' at time t'. The functions -D+(-) 
are initially right (left)-propagating 5 pulses originating 
at point x' at time t' . 

Let us for simplicity assume that the two supercon- 
ducting contacts are pointlike in the low-energy model, 
that g — 1, and that the helical liquid is homogeneous. 
Then, D±(x, t; x', t') = 9(t - t')S(x - x' T v F (t - t')) and 
we find 



i(x,t)= [e(x- Xj )i T {i +tj (t- 

i=i,2 



X 



6( Xj - x)l T (I-,j(t 



v F 
\ x ~ x j 

VF 



')) (38) 



))]• 



i 9 



I-j = 2iT 



+ r + _v+'0- + h-c \ x =xj 

2ir ++ f/; + V+ + «r + _?/; + ?/)_ + h.c. 
ip-ip- + ir^ ip+ip— + h.c. 



(39) 



(40) 



where It is the small contact size, and the expectation 
values (•) closely correspond to the different parts of the 
injection currents Ij y si/S2 evaluated in the previous sec- 
tion. Indeed, 



tf-,i(*)>=ijf(*) + ^sf(*)» 

(t,2(i)}=^ + (t) + ^ 2 -(*), 
(I-,2(t))=js2(t) + ^j+z(t). 



(41a) 
(41b) 
(41c) 
(41d) 



The physical interpretation is particularly simple: the 
contacts at x\ and X2 inject current to the helical liq- 
uid. The component due to H — tunneling splits evenly 
to the left and right-moving modes, whereas the ++ and 
components end up solely in the + and — modes, re- 
spectively. Within each edge mode, the injected current 
propagates with the Fermi velocity, as indicated by the 
retarded time arguments. 



The calculations done in the previous section indicated 
that in this case jg-[~ = and jg^ = to leading order in 
a. Therefore, essentially all of the current injected by the 

++ and tunneling in fact flows only to the reservoirs 

that maintain the chemical potential of the helical liquid 
at fi = 0, rather than between the two superconducting 
contacts, which can be verified by computing the current 
at x < x\,X2 and at x > x\,X2- The effect essentially 
amounts to a modulation of the NS conductance between 
the superconductors and the normal leads by the (time- 
dependent) phase difference <pi(t) — if2(t) between the 
superconducting contacts. 

Based on the above results, we can write down an ex- 
pression for the part of the NS current [see Fig. |ljc)] 
that depends on the phase difference, in the configura- 
tion Vi = V 2 = V: 



81ns = 



2l T \X ++ \ 
3tt 



Vcos{2Vd + 4> ) cos(ip ) (42) 



x [(V/A) 2 + 4tt 2 (T/A) 2 



Al T \X+. 



sin(2F<2) cos(^o) — 



smh(2zd) 



Note that the modulation of the NS conductance from 
the H — channel decays exponentially as the tempera- 
ture increases, whereas the ++ contribution does not. 
The same situation should persist in all orders of pertur- 
bation in the effective Hamiltonian for the H — tunnel- 
ing: the terms coupling to ipo contain inequal numbers of 
ip+(d)ip-(d) and ?/>t(d)?//L(d), which implies that the cor- 
relation function is of the form [B+(d,t)B-(d,t)] x 0(1) 
and thus has an overall exponential prefactor e ~ 2 ^ Td / v F _ 
Therefore, there in principle is a temperature regime at 
T 3> hvp/kgd in which the leading contribution to the 
tpo dependence of the NS current comes mainly from the 
++ tunneling, despite the power-law suppression of this 
channel in helical liquids. The physical reason for the 
difference can be seen in Fig. [2j for the H — channel 
an electron pair injected to energies ±e and traversing 
through the junction obtains an energy dependent phase 
factor e l ( k i+ k 2)d = e -2ied/v F ^ typica i f Andreev reflec- 
tion, which averages towards zero when a finite energy 
window is considered. For the ++ channel, because of the 
linear spectrum, the corresponding phase factor e 2lkpd is 
energy-independent and no such averaging occurs. 

The above result requires validity of the perturbation 
theory, i.e., It^-\ ~ Rk/^Rn ^ 1 where l T is the con- 
tact length. If this condition is not satisfied, an addi- 
tional contribution decaying only as 1/T in temperature 
arises in the H — channel!^ This additional proximity 
effect contribution is similar to what occurs in metallic 
systems of a more macroscopic sizepSl although there it 
can be much amplified as the electrons can stay a long 
time near the NS interface due to impurity scattering. 

With finite repulsive interactions (g < 1), also the H — 
contribution to the NS conductance obtains a power-law 
prefactor according to the scaling dimensions, 1^ oc 
(E/A) 2 /b- 2 with E = max(T,V), and the prefactor of 
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the ++ part is modified, /++/— oc {E / Af^+9 x )-2. 
According to the correlation functions (31), exponential 
decay will also appear in the ++ part due to charge frac- 
tionalization, but it will be weaker than in the -\ — part 
for all values of g. 

A second distinguishing feature of the ++ contribution 
to the NS current is that it is expected to oscillate not 
only as a function of the bias, but also as a function of 
the Fermi wave vector appearing in the dynamical phase 
(j)Q. In metals or other systems where kp is large, the 
wavelength of such oscillations would be on the atomic 
length scales, and the contribution would average to zero 
[as ~ sinc(fc_F«;) 2 ] over any practical contact size w^l 
However, this needs not be the case in HgTe-QW (or in 
nanotubes, see Ref. [TO)) when the Fermi level lies close 
to the Dirac point: for example assuming \fi — -Ebiracl ~ 
M/6 ~ 1.5 meV one finds 1/kp ~ 150 nm. Such length 
scales are likely experimentally accessible. 

One should also note that the finite wave velocity com- 
bined with the ac Josephson effect causes some additional 
effects. The current propagates at the (renormalized) 
Fermi velocity vp/g, rather than at the substantially 
higher speed of light c at which electromagnetic exci- 
tations propagate. Assuming only the H — channel con- 
tributes, one can find the spatial dependence of current 
between the two contacts: 

I(t) cx sm[2Vt - — ] + sm[2Vt - 2 Xi^A ] . (43) 



v F 



Vp 



Based on this, it is clear that for biases V > hvp/ed be- 
tween the two superconducting electrodes, the ac Joseph- 
son effect must be associated with appreciable standing 
wave oscillations in the charge density. This behavior is 
not specific to helical liquids: a similar spatially resolved 
calculation as above for the spinful liquid ac Josephson ef- 
fect of Ref . US] should also produce this feature. Whether 
such effects are observable in reality, however, depends 
on how realistic the model assumptions about screening 
are in the systems studied (see also Ref. [31]) . 



V. DISCUSSION AND CONCLUSIONS 

In this work we considered the proximity effect induced 
in a helical edge state, taking into account a spatially and 
energetically non-constant spin quantization axis. Such 
rotation of the spin axis naturally arises from the spin- 
orbit interaction in real materials such as the HgTe QWs, 
for example in a controlled way by structure inversion 
symmetry breaking Rashba terms. This has the conse- 
quence that the singlet correlations in an s-wave super- 
conductor can also induce a proximity effect in the same 

channel of left and right-movers (r ++ / ) in addition 

to the usual term where the correlation is between op- 
posite chiral states (with amplitude ). We derived 

a description of the proximity effect in both channels in 
the presence of Rashba interaction using a simple model 
for the tunneling between the superconductor and the 



helical edge state, respecting spin conservation and time 
reversal symmetries. 

The extra transport channels (+ + / ) describe pro- 
cesses that are in principle parasitic for the splitting of 
a Cooper pair into two electrons propagating into differ- 
ent directions (the H — channel}^. For a single super- 
conducting contact to the helical liquid, the scaling with 
temperature (or bias voltage) at low energies however al- 
ways favors the H — channel. In Ref. [32 the two-particle 
tunneling into the bulk of a spinful Luttinger liquid was 
found to be suppressed in a power law in 1/ A similarly as 

here, but there the tunneling into the + + / channel 

was found to be dominant. The difference arises because 
in a spinful liquid the two opposite spins can tunnel into 
different spin channels, and therefore no Pauli-blocking 
factors appear. 

Observing effects related to the same-mode tunneling 
(r ++ ) is likely rather challenging, as they can be sup- 
pressed relative to r_| by several factors: the power law 

suppression (T/A) 2 + (V/A) 2 from exclusion principle, 
suppression of the tunneling factor T ++ itself, and aver- 
aging effects related to contacts if they are larger than 
1/kp (i.e. ~ 150 nm for parameters in Fig. [5]). How- 
ever, by observing the dependence of the NS conductance 
on the superconducting phase, the relative difference can 
be reduced due to the exponential dephasing of the H — 
contribution at high temperatures. (In the case that the 
only mode of transport is via the H — channel, the mod- 
ulation would still contain features distinct to ballistic 
transport, such as oscillations as the bias voltage is in- 
creased.) The question is therefore more on how small 
signals can be detected in the conductance, oscillating 
with the phase difference ipo, and how large the thermal 
factor 2^k B Td/hvp can be made before inelastic interac- 
tion effects (e.g. electron-phonon scattering), which we 
have neglected, start to play a role. 

We find that controlling the spin axis via external elec- 
tric fields in HgTe-QW in general requires very strong 
fields, because of the weak coupling of the additional 
spin-orbit interactions to the edge states. For our case, 
this makes achieving a large T ++ more difficult, and 
may in general pose problems to proposals relying on 
the control of the spin axis. Making an optimistic es- 



timate, we find from Eqs. (23) and (42) that the ra- 



tio of the two contributions to the amplitude of phase- 
dependent oscillations in the conductance is (V <C T, 
g=l,E z ~ 100 m v/ nm ) 



5G ++ ((p Q ) 



27T 2 

IT 



k B T 
A 



sinhc 



2Trk B Td 

hvp 



(44) 



k B T 
M 



sinhc(27r ksTd/hvp) . 



where sinhc(a;) = smh(x)/x. With finite e-e inter- 
actions [cf. Eq. ( |B25 1], the ratio is multiplied by 
(A/T) 2 ~ 2s , making the result depend only on A/M in 
the limit g — > 0, and the exponential dependence be- 
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comes ~ exp(2irg[2 — g]kBTd/Hvjr). Taking junction 
length d = 3 /jm, the temperature scale of the exponen- 
tial suppression factor is = hvp/(2Trd) s» 0.2 K, and 
the ratio becomes unity at the cross-over temperature 
T, ps E T {\og[2(M/E T ) 2 } - loglog[2(M/S T ) 2 ]} « 2.2 K, 
which depends weakly on M (here M = — lOmeV). 
Given a suitable superconducting material, this should 
be achievable. 

Another option for amplifying the same-mode tunnel- 
ing could be to break the time-reversal symmetry and 
introduce additional spin flips or spin rotation, for exam- 
ple via magnetic impurities or ferromagnets. The effect 
could still be detected in the NS conductance, as that 
conclusion is only based on the generic form of the low- 
energy effective Hamiltonian. 

Observe that in our analysis the true 2D nature of 
the edge states in HgTe-based QWs was important. The 
( + + / ) proximity channel cannot be found in a com- 
pletely ID description, as in such a picture the spin quan- 
tization axis is simply rotated globally by the Rashba 
terms (cf. Refs. fl3l and l33|). Such rotations can have 
no consequences for Cooper pair injection into a single 
edge, due to the s-wave symmetry of the pairing [cf. 
Eq. jjjh]. Spatially inhomogeneous Rashba interactional 
could, however, induce a finite + + / amplitude. 

As in other systems with small critical currents,^ 
also here thermal fluctuations in the superconducting 
phase difference are a problem for measurements of the 
temperature-dependence of the Josephson effect: the 
temperature scale relevant for the phase fluctuations, 
Ej — hl c /2e, is smaller than the intrinsic one, Ej> — 
hvpjd. More complicated measurement schemes^^SI 
than the simple current-biased setup in Fig. [ljb) may 
nevertheless help in overcoming this problem. One 



should, however, note that only the Josephson current 
is a problematic observable in this respect. The mea- 
surement of phase oscillations of the NS conductance in 
the setup of Fig. [ljc) is expected to suffer much less from 
phase fluctuations, as there the phase difference is locked 
by the magnetic flux and the large critical current of the 
superconducting loop itself. 

In summary, starting from a tunneling Hamiltonian, 
we derived an effective low-energy theory describing the 
superconducting proximity effect in the helical edge state 
of a 2D topological insulator. We showed that in these 
systems, despite the s-wave symmetry of the supercon- 
ductor, correlations can occur both in (+ + / ) and be- 
tween (H — ) the left and right moving modes, and within 
a simple model, we estimated the expected magnitudes 
for the effective proximity gap parameters in HgTe/CdTe 
quantum wells. Based on the effective Hamiltonian, we 
studied the dc and ac Josephson effects in the helical liq- 
uid, and considered phase-dependent oscillations of the 
NS conductance. In nonequilibrium, we found that cor- 
relations within the same mode can give rise to a long- 
ranged interference effect, which could act as a signature 
of their presence. Our results also shed light on the mean- 
ing of " spin" in the helicity of these edge states which is 
of importance if one intends to use these edge states for 
spin-injection or spin-detection. 
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Appendix A: HgTe/CdTe QW edge states 

The edge states of a HgTe-QW can be described within 
the four-band model introduced in Ref.;2, Here, we derive 
explicit analytical expressions for the edge states in a 
single edge following the approach of Ref. |2ll for use 
in Section. |IIIB[ and to demonstrate that the direction 
where the 4-band spinors point is independent of k and 
M, in the absence of inversion symmetry breaking terms. 

The four-band Hamiltonian reads 



_ (h(k) 

a ~ \ o h{-ky 



(Al) 



h(k) = e(k)a + d(k) ■ a , e(k) = C- D\k\ 2 , (A2) 
d(k) = (Ak, -Ak y , M - Bk 2 ) . (A3) 

For the parameters A, B, C, D we use values from 
Ref. HQ A = 365meVnm, B = -706meVnm 2 , D = 




(b) 



r 



r 



FIG. 7. (a) Cooperon giving a contribution to the Josephson 
current, (b) Integrating out the relative coordinates gives an 
effective description of Andreev reflection. 



— 532meVnm 2 and take C = (it only shifts the Dirac 
point). For an edge with QW lying at y > 0, with the 
wave function vanishing at y = 0, the edge eigenstates 
are: 







(A4) 



$ Q , fe - Ne~ ikx [e- x ^ - e- x ^] ( ) , (A5) 



-y/D-B 



where 



-DM , AVB 2 - D 2 
E = — ak- 



B 



B 



Ai/2 = V fc2 + f t Vf 2 -q 2 , 



(A6) 
(A7) 



F= A2 Z 2 i BM ± DE \ Q > = ^r, (A8) 



2{B 2 - D 2 ) 



B 2 - D 2 



and N is a normalization constant. 

Note that the spinor points to a single direction inde- 
pendent of k or energy E, so that the above results are 
of the form 



* +,k = ( n ) i+,k > = ( I ) f-,k , (A9) 



where the constant spinor \ i s normalized (x^X = 1) an d 
depends only on the parameters B and D. The envelope 
f a ,k(x,y) = N'e lkx [e- Xi y - e- X2V ] is a scalar function. 
Using the above parameters we find x = (0.35, 0.94) T , 
so that the spinor has the main contribution in the HI 
band. 



Appendix B: Low-energy Hamiltonian 

In this Appendix, we derive the effective Hamilto- 
nian in Eq. Q via perturbative renormalization group 
(RG)p3 aiming to approximate the Cooperon term [see 
Fig.[7]ja)] appearing in the perturbation expansion in the 
tunneling with the effective term in Fig.[7]jb). When com- 
puting the Josephson current, for the channel this 

approach is compatible with that used in Ref. [18] and 
elsewhere in the long-junction case. For the T ++ and 

r channels we however need to pay more attention to 

the tunneling elements. 
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We take our effective Hamiltonian to have the form: 
H cB = H Q + H T + H T2 = H Q + V, (Bl) 
//■/ > / da; d 3 r't aa ,(x, ?)il>l t (x)ip s *> (»*) + h.c. 

(B2) 



Ht2 = 



a/9 



dec r a p(x)tj) a (x)ipp(x + a) + h.' 



(B3) 



and rescale only the cutoff in the helical liquid in the 
progress of RG. Because correlations in the superconduc- 
tor decay exponentially at distances > A -1 , and because 
the superconducting gap prohibits dissipation at low en- 
ergies inside the superconductor, Ht can be neglected in 
calculations after the scaling to low-energy length scales 
^> A -1 is done — which reduces the effective Hamilto- 
nian to that in Eq. Q. 
The scaling equations read 



dfc. 



dl 
dl 



[2 - m]t aa , (B4) 
[2 - V2,a0]r a p + S a p(l) , (B5) 



where r\\ = (g + g~ r )/A is the scaling dimension of 
e i<P± appearing in tunneling Ht, and f]2, a -a — l/<7 and 
V2,a a = 9+9 1 are scaling dimensions of the operators in 
Ht2- Although essentially a standard calculation, below 
we explain the derivation of S a p in detail. 
Below, we need the factorization 

i> ai (xi,Ti)i> a2 (x2,T 2 ) = (a q Q ) 2 ^ UaiiT ' )Ua2iT2) (B6) 

Z7TO0 

v ik F (a 1 x 1 +a 2 x 2 ) 



X2,n - t 2 ) 



where q$ — 2ir/L is the infrared cutoff, and the correla- 
tion functions read 

C ++ (z) = ( qoZ r 1 ^ 2 ( qo z*y- 1 ^ 2 (B7) 

C a ^ a (z) = \qoz\ {1 ' 92)/{29 \ (B8) 

where C a p(x,T) = C a p(z), z — vft — ix, and C (z) = 

C ++ (z)*. Observe that C aiOL2 {z) cx qQ 2,aia2 ~ Vl . 

To perform the RG steps, we also need the correspond- 
ing operator product expansions. Taking sign changes 
due to Klein factors and time ordering into account, we 
find [cf. Eq. pi])]: 



T[ip ai {zi)ip a2 (z2)] 

l s , f zi + z 2 ,z 1 +z 2 

u aia2 {zi - z 2 )ip ai { — - — )rp a3 { — - — 



where 



!> i(a 1 -a 2 )k F x/2 



sgn(r)^ 



(B9) 

oo) , 

(BIO) 



C QljQ2 (xsgn(T), |r|) 



C a 



2 (0,o ) 



and S ai a2 in the sign factor arises from the fact that 
U ai U a . 2 = (-l) 1+s °i-°>U a2 U ai . 

The source term S a p(l) for the Andreev reflection pro- 
cesses T a p appears from the second-order term in the 
pertubation expansion of the partition function, Z/Z = 

(Te-f° dTXV ^)o = l + c 1 \ + c 2 \ 2 + .... Combining two 
Ht and using the operator product expansions gives a 
contribution to T a p. We also trace out the supercon- 
ductors at this step, factorizing the expectation value to 
(. . .) = (■ • -)hll,o(- ■ -)s,o- This yields the result 

dc f 

— j- = / d 2 z ^2(T[i/> a (z)i/>p(z + a )]) S a /3(l,z), 

(Bll) 
(B12) 

(B13) 
(B14) 



S a p(l,z) = a [d r f a p(l,z,r)] r = ao , 
fap(l,z,r)= / d 2 z'c a p(l,z, z') , 

J\z'\<r 



C Q /3 



(M,z')^/dMdM E 



CriCT2=t4 
/ /o £>\*4 ( ™ ™' /I ^ \* 



x t aai (x + x'/2, fx)*t^ 2 (x - x'/2, f 2 ) 
x F\cr 1 ,r l 1 ,T';a 2 ,r' 2 ,Q)u a p{l,z') , 

where d 2 z is shorthand for dxdr . This closes the set of 
equations. 

We can now solve the scaling equations: 



tester' (0 ^ ^ ^ ^a<r(0) ) 



(B15) 



T aP {l,z)= d s ^ 2 -^'^ l -^S a p( Sl e l ' s z). (B16) 



The integral appearing in T a p can be simplified by sub- 
stituting in the scaling obtained for t a p , and undoing the 
rescaling of length scales in the remaining integrals. This 
yields: 



u a p(l,z] 



— J2vi-r) 2 , a f)]l 



^u aP (0,e l z) 



(B17) 

c a0 (l, z, z') = e^-^ l c a0 (O, e l z, eV) (B18) 
f a p(l, z, r) = eP-fr.-d'/a/KO. e l z, e l r) , (B19) 
S aP (l,z) = e^-^^ l a [d r f a p(0,e l z,r)} r=aoe i (B20) 

And further, 

T al3 (l, z) = el 2 -*-*]' f d 2 z' c aP (0, e l z, z') , 

J a <\z'\<a(l) 

(B21) 

where a(l) — a Q e l . Here, c a p decays fast for \z'\ > |A| _1 
due to the decaying F functions and the assumedly short 
range of tunneling. Therefore, at long length scales 
a = aoe 1 3> |A| _1 we can replace the upper limit in 
the integral: a t— > oo. 

Undoing all length rescaling, we can write the result in 
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the form of Eq. Q , with 



Within these assumptions, we get 



T aP {x) = * / d 3 r[d 3 r' 2 dx'dT'e l{a -^ kFX ' ? 2 (B22) 

x^+^i^-y,^ fa - r'; f a , 0) 
a l -^C a p{x' sgn(r'), |r'|) sgn(r')^ 



l-»7a,ai8 



C Q) g(0,ao) 



where 



Consider the noninteracting case, .9=1. There, 



,l-r/2,c3 



C«^(a;sgn(r), |r|) sgn(r) 6 



a [vft — ix\ 



C a p(0,a o ) 



(B24) 



respectively. One also notes 



+ 



so we re- 



fer a/3 = H — h and a/3 = 

that r = — r_i 

define i— > 2T^ as the sum of the two and drop 

the — h term. Finally, going into Fourier representation 
yields Eqs. ^ and §6§. Due to the integrals over x' and 
t' extending over the whole range, and the correlation 
functions C a p being local in frequency and energy, only 
certain energies and momenta contribute in the final re- 
sult. 

To find out the effect of interactions, one needs to 
roughly estimate the result from Eq. (B22). First, 
since fc^g in the superconductor is a short length scale, 
we take F(r[, r'; f 2 , 0) i-> ^(t')^ - rj), F{t') = 
iV>Aif (|r|A), where N F is the normal-state DOS at 
Fermi energy in the superconductor, and \t\T <C 1, 
T <C A <C \M\. We consider tunneling that is local on 
length scales of 1/A (and 1/T and the other low-energy 
scales) and replace P^ ~ JCS(x')S(x — x[) with JC a con- 
stant. Based on the model in Sec. ( IIIB[ ), the Fourier 

transform of P++ in x' satisfies P++ ~ —ik(z$/M)P^ 

on long wavelengths \k\ <C M, with P + _ constant in k. 
In real space we then have P++ ~ JCjjd X 'S(x')5(x — x'i). 



KNpa^a^A I dr' ir (A|r'|)C+_( — ) 



ICN F {aA) 1 - 1 ^(a A) s± ^- 1 q{ 



9 - 9 s 



P ai a 2 (xi,r' 1 ;x 2 ,r 2 ) = [t ctll (x ll r f 1 )t a2t (x 2 ,r 2 ) (B23) 
- t ai ^{x u r l )t a2i {x 2 , r 2 )} + [r[ <H> r 2 ] , 



and we have made use of the singlet symmetry of the F 
function. The cutoff a in the theory specified by Eqs. Q 
and (B22| can be chosen freely, but taking a — |A| 
is natural as the source term in the original RG stops 
contributing at that length scale. 



V F T - IX 

x Ida,- ImC ++ ( )\x'=o 

9o 

!^(aA) 1 -»- 1 /» (fl0 A ) <-i_^ 



9o 
(B25a) 



(B25b) 
(B25c) 



v F \M\ 



(B25d) 



x 9(^-D 



where q{x) = 2 X - 1 T(±4^) 2 ~ 1. The effective tun- 
nel rates obtain an identical scaling in the bare short- 
distance cutoff a Q related to interactions, which appears 
because of a renormalization of the tunneling elements 
t aa . The low-energy scaling with a follows the scaling 
dimensions in the effective Hamiltonian. Finally, T++ 
has an additional factor ZqA/vf\M\ that is a signature 
of the Rashba coupling. 

The expression ([5| for T ++ deserves some comments: 



First, we know that F(lj) 



1 



for 



0, so 



that <9 w P| w= o vanishes, and from Eq. (11) we know that 
P++(—k) = — P ++ (fc) which means that dkP++ is even 
in k and can be finite at k = 0. Note that the gradi- 
ents , dk appear because the boson correlation function 



[see Eq. (B24)], 



( e #+( a .- r )e^+( ' ))o vanishes at x,r 
which reflects the fermionic exclusion principle. This is 
the reason why the effective Hamiltonian contains a term 
resembling more tp+^Pi- than tp+kip+, which is in agree- 
ment with the results of Ref. HTH 

Finally, we observe that in the noninteracting case, 

is related to the leading-order off-diagonal Nambu 

component of the self-energy. The factors T ++ (and 

in the interacting case) however in general contain 

additional information, as the out-integration of short 
length scales captures the renormalization from interac- 
tions, and the effect of the exclusion principle when av- 
eraging ip(x + x')tp(x) over short distances x' . 



Appendix C: Current operator 

For completeness, we include here a derivation of 
Eq. (36 1 that shows the result obtained in Ref. [55] ap- 



plies also to time-dependent perturbations. Related re- 
sults can be found e.g. in Ref. [37l and a special case of 
the present result is given in terms of path integrals in 
Ref.HSJ 



14 



Consider the Heisenberg equation of motion under a 
Hamiltonian H = Ho+V(t), where Hq is given in Eq. ([lj, 
and the perturbation V(t) is switched on at t > 0. Iter- 
ating the equation of motion for d x <f> twice, one obtains 

d 2 t (d x <p) - d x (u 2 d x (d x <t>)) = s(x, t) (CI) 

s(x, t) = [H , [H , d x ct>}} - [H, [H, d x 4>]] + i[H, d x <f>] , 

(C2) 

where H = V contains the explicit time dependence of 
the Hamiltonian, and u{x) — vp/g(x) is the renormalized 
wave velocity. The solution to this linear equation can 
be written in terms of the retarded Green function of the 
wave equation on the LHS: 



d x (j>(x,t) = d x (j) (x,t) 



(C3) 



Ax' C R {x, t; x', 0)i[V, d x <t> (x' ', 0)] 
<1/' / dx'C R (x,t;x',t')s(x',t'), 



where d x (f>o evolves under Hq, and the second term en- 
sures that the initial condition dt(d x (j>) = i[H,d x (j>] is 
satisfied — this follows from d t C R (x, t; x' , i')| t _ >t / +0 + = 
S(x-x'). 

We can also rewrite s using properties of the fields and 
H : 



s(x,t) = -d x (v F g(x) 



_ 2 SV(t), 



5(j){x)> 

where we noted the correspondence 
SA 



sv(t) 

5$(x) 



[A,a x 4>[x)\ 



5d{x) 



[A, d x $(x)} 



SA 



(C4) 



(C5) 



valid for functionals A — A[d, <fi]. The expression for s 
can be substituted in Eq. ( C3 1 , and an integration by 
parts transfers the gradients to operate on C R . One of 



the resulting boundary terms cancels the second term in 



Eq. (C3), and the others vanish, provided the perturba- 
tion s vanishes at x 



±oo. 

We then find the exact result 

poo 

d x (f>(x,t) = d x ef> eq (x,t) + J <!/' / <l.r' 



(C6) 



v F g(x')- 2 id x ,C H (x,t;x',t')] 



5<f>(x') 



[d t ,C R {x,t;x',t')) 



8§{x') 



We can simplify this further by making use of prop- 
erties of the 1-D wave equation. The Green function 
C R (x,t;x',t') satisfies the initial value problem 

[d 2 t - d x u 2 d x ]C R = , at t > t' , (C7) 

C R = , d t C R = S(x -x'), at t = t'. (C8) 
If u is a constant, the solution is a sum of two wavefronts 
C R = C R +C R , C R = {Au)- 1 9{t~t')sgn[±(x'-x)+u(t- 
t 1 )]. This is valid in the limit t — > t' also if u is smoothly 
spatially varying — the wave equation only sees u around 
x' . Since the wave equation is linear and its solution is 
unique, the Green function can always be decomposed 
to these two parts. Let us now define D± = 2dtC R and 
F± = 2d X 'C±. They satisfy the wave equation at t > tf, 
and the initial conditions are inherited from the t — > t' 
behavior of C±: 

D± = 8(x - x') , d t D ± = tu{x')S'(x - x') , (C9) 
F ± = ±u{x')- 1 5{x - x') , d t F ± = -S'(x - x') . (CIO) 

Due to linearity, clearly F± = ±u(x')D±. Be- 
cause C R (x, t; x' , t') — C R (x,x',t — t'), we then find 
d t ,C R = -6 t C R - -^±±^ and d x ,C R = ^±+^ = 
u(x')^ 



ing j± 



Substituting these to Eq. ( C6 ) and defin- 



y Sep 



± 



SVi — sv 



81) - 



, we arrive at Eq. ( 36 1 



